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A long standing problem of domain switching process - how domains nucleate - is examined in 
ultrathin ferroelectric films. We demonstrate that the large depolarization fields in ultrathin films 
could significantly lower the nucleation energy barrier (£/*) to a level comparable to thermal energy 
(fesT 1 ), resulting in power-law like polarization decay behaviors. The "Landauer's paradox": U* is 
thermally insurmountable is not a critical issue in the polarization switching of ultrathin ferroelectric 
films. We empirically find a universal relation between the polarization decay behavior and (7*/fcsT. 
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Recent advances in complex oxide thin film synthe- 
sis and in first-principles calculations have intensified 
the basic research on ferroelectricity at nanoscale di- 
mension 0, Q • The cooperative nature of ferroelectric- 
ity is expected to induce different polarization states 
at nanoscale due to low-dimensional boundary condi- 
tions and the strong interactions of polarization with 
strain,charge, and other electromechanical parameters 
IH S 0, El Many workers have reported in- 
triguing physical phenomena occurring in ultrathin fer- 
roelectric (FE) films and other nanostructures, such as 
intrinsic size effects y, |3J , strain-enhanced FE proper- 
ties Q , unusual low-dimensional phases [f| , and domain 
patterns |g, |j| • On the other hand, the mechanism and 
domain dynamics of FE switching in ultrathin films have 
rarely been investigated in spite of their scientific and 
technological importance H |tj El • 

Historically, the mechanism of polarization switching 
dynamics in ferroelectrics has been the subject of a great 
deal of research. It is now believed that the polariza- 
tion switching takes place not by coherent switching, but 
by the nucleation and growth of new domains 0, 0, El • 
However, there still remains an unsolved issue: How do 
the domains nucleate? In the late 1950s, Landauer em- 
phasized that a thermodynamic nucleation process can- 
not play a role in FE domain switching [TlT ] . For a nucleus 
with reversed polarization to be created by a thermody- 
namic process, an energy barrier for nucleation U* should 
be thermally overcome, as shown in Fig. 1(a). However, 
Landauer's and later estimates showed that U* is practi- 
cally insurmountable by the thermal activation process: 
U* > 10 s k B T at an electric field E ~ 1 kV/cm (a typi- 
cal value of the coercive field for bulk ferroelectrics) [llj 
and U ~ 10 3 k B T at E ~ 100 kV/cm (a typical value of 
the coercive field for most FE thin films) floL fl3| , where 
ks is Plank's constant and T is temperature. According 
to these estimates, it is difficult to understand the ob- 
served domain nucleation in ferroelectrics. This problem 




FIG. 1: (color online), (a) A schematic diagram of electro- 
static energy for the nucleus formation ([/) as a function of 
the nuclear radius (r). Schematic diagrams of (b) the half- 
prolate spheroidal nucleus formation and (c) the cylindrical 
nucleus formation, with reversed polarization. 



has been known as the "Landauer's paradox". To over- 
come this difficulty, numerous workers have assumed that 
the nuclei could be formed inhomogeneously due to ex- 
ternal effects, such as defects 

.12]: long-range interaction 
between nuclei |l3j and FE-electrode coupling 14] . How- 
ever, in this Letter, we will demonstrate that U* could 
be thermally overcome in the ultrathin FE film. The 
thermodynamic nucleation of reversed domains due to a 
large depolarization field Ed can govern the polarization 
switching process, without inclusion of any external ef- 
fects, in ultrathin epitaxial BaTiOs films. 

Fully strained SrRuOs/BaTiOs/SrRuOs capacitors of 
high quality were fabricated using pulsed laser deposition 
with reflection high-energy electron diffraction jl]. The 
thickness of BaTiOs (rfsro) is between 5.0 and 30 nm. 
We were able to directly measure P-E hysteresis loops 
and polarization decay behaviors even for a 5.0 nm thick 
BaTiOs film [J, [lj| UJj • We focused on the polarization 
decay behaviors for different BaTiOs film thickness and 
measurement conditions. As shown in Fig. 2(a), these 
SrRu0 3 /BaTi0 3 /SrRu0 3 capacitors with a 10x10 /j,m 2 
area displayed an intriguing time-dependent polarization 
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FIG. 2: (color online), (a) Time-dependent net polarization 
changes, AP(t), without any external field (E ext ). (b) AP(i) 
of the 5.0 nm thick BaTi03 capacitor under numerous val- 
ues of E £X t- The left- and right-hand diagrams of the inset 
show electric pulse trains used to measure non-switching and 
switching polarizations, respectively. AP(t) corresponds to 
the difference between these polarizations, (c) The relation 
between n and the total internal electric field (namely, |-Ed|- 
l-Eeitl)- The open and solid symbols correspond to the data 
with and without E ext , respectively. 



decay behavior. We measured the time-dependence of 
polarization decay using electric write /read pulses sep- 
arated by a given decay time (t) 0, ^| : the right and 
the left schematics in the inset of Fig. 2(b) show the 
electric pulses used to measure the switching and non- 
switching polarizations, respectively. Figure 2(a) shows 
that the polarization decay follows a power-law behavior 
at room temperature [l^: namely, AP(t) ~ t ~ n . As 
the film becomes thinner, the value of n increases sub- 
stantially. In particular, for the 5.0 nm thick capacitor, 
AP decays below 1% of its initial value within 10~ 3 s. 
This rapid polarization decay must be caused by domain 
backswitching upon removal of the electric field and is 



the subject of the discussions that follow. 

Within the capacitor geometry, the discontinuity of 
polarization at the FE/electrode interfaces induces po- 
larization charges. The free carriers inside the electrodes 
compensate the resulting charge deficit. However, the fi- 
nite value of the Thomas-Fermi screening length for the 
free carriers makes the compensation incomplete, result- 
ing in Ed inside the FE layer ^t|- It should be noted 
that the generation of Ed is inevitable, which will limit 
the performance of most FE capacitor-type devices. 

To obtain further information on the detailed effects 
of Ed on the polarization decay, we measured the time- 
dependence of AP with an external static field E ext , ap- 
plied in the opposite direction to that of Ed during the 
electric pulse measurements, displayed schematically by 
the red dashed lines in the inset of Fig. 2(b) [l^. Since 
E ext will partially cancel out Ed, it will slow down the 
polarization decay and consequently reduce the values of 
n. As shown in Fig. 2(b), n decreases significantly with 
an increase of E ext and should finally become zero when 
E ext becomes nearly the same as Ed- This E ext value at 
n ~ provides an experimental value for Ed |15j . For 
the 5.0 nm thick BaTiC>3 film, Ed can be as large as 800 
kV/cm. Note that this value of Ed is even larger than 
the coercive field value, i.e. about 300 kV/cm, which was 
measured from direct P-E hysteresis measurements 

Figure 2(c) shows the dependence of n on the total 
internal electric field, i.e. \Ed\-\E ext \, for the BaTiOa ca- 
pacitors. The abscissa-values of the solid symbols, which 
were measured with E ext —0, correspond to the experi- 
mentally estimated Ed values. They clearly show that Ed 
becomes larger as dsro becomes smaller. As the \Ed\- 
\E ext \ increases, the values of n increases for all samples. 
However, we could not find any simple universal relation 
between n and |2£d|-|.E ea ; t |. 

Further insights could be obtained by looking into the 
mechanism of domain switching dynamics in an ultrathin 
FE film. Note that decrease of AP means an increase in 
the number of FE domains with a reversed polarization. 
It is generally believed that the growth of the opposite 
domains occurs in three processes, namely (a) formation 
of nuclei with opposite polarization, (b) their forward 
growth, and (c) sidewise growth (called, the domain wall 
motion). For epitaxial FE films thicker than 100 nm, the 
domain wall motion has been reported to be very impor- 
tant [l^. However, for ultrathin FE films, domain wall 
motion becomes extremely slow: for example, the speed 
of the domain wall motion was estimated to be about 1 
nm/s for a 29 nm thick Pb(Zr,Ti)0 3 film 0. In addition, 
the forward growth of the nuclei is known to occur very 
fast (i.e. a speed of 1000 m/ s) 9] , so the very slow domain 
wall motion could not play a major role in the polariza- 
tion decay behavior observed in Fig. 2(a). Therefore, 
we propose that, in the ultrathin FE films, the rapid po- 
larization decay should be caused by nucleation-governed 
domain dynamics and that the large Ed value could make 
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FIG. 3: (color online), (a) Normalized AP from the exper- 
imental data (the solid symbols) and the Monte-Carlo sim- 
ulation (the solid lines). The black dotted line in the inset 
shows the Monte-Carlo simulation results for the 5.0 nm thick 
BaTiOs capacitor without lateral screening of the interaction 
term between long-range domains in the Hamiltonian. (b) 
Normalized AP of the 5.0 nm thick BaTi03 capacitor un- 
der numerous values of E ex t- The solid symbols and lines 
correspond to the experimental data and the Monte-Carlo 
simulation results, respectively. 



the thermodynamic nucleation process possible. 

Figures 1(b) and 1(c) schematically display forma- 
tion of nuclei with half-prolate spheroidal and cylindrical 
shapes, respectively, under an electric field. Here, d is the 
FE film thickness, and r* and I* represent critical radius 
and length of the half- prolate spheroidal nucleus, respec- 
tively. Kay and Dunn calculated the electrostatic energy 
U for such nuclei formation [l^, and showed that there 
should exist the nucleation energy barrier U*, as shown 
in Fig. 1(a). When I* < d, the spheroidal nuclei are likely 
to be formed, and vice versa. They also showed that U* 
is proportional to a w 3 /E 5 / 2 for the half-prolate nucleus 
and a w 2 /E for the cylindrical nucleus |l8| . where a w is 
the domain wall energy. Using a w = 10 mJ/m 2 0, the 
crossover from the half-prolate spheroidal to cylindrical 
shapes will occur at d ~ 10 nm |2fJ. Note that the very 
large Ed values of our ultrathin BaTiOa films will make 
U* small enough to allow a thermally-activated nucle- 
ation process. The experimental Ed values of our BaTiC>3 
films are between 300 and 800 kV/cm 0. With these 
Ed values, the U* /ksT values are estimated to be be- 
tween 4 and 20 at room temperatures, which are easily 
accessible by thermal energy. 

To explain the observed power-law decay behavior, 



we used numerical simulations adopting the Monte- 
Carlo algorithm, which takes into account the thermally- 
activated polarization reversal process. We approxi- 
mated the film to be composed of 128 x 128 single do- 
main cells lying on the a6-plane, whose size is 2.0 nm and 
height is the same as the film thickness, with a periodic 
boundary condition. Each cell has a uniform P, either 
pointing up or down along the c-axis. We used a dimen- 
sionless parameter x (= P/P s ), and set the value of x to 
be either +1 or -1 for the bi-stable states, where P s is 
spontaneous polarization. During the switching process, 
the Hamiltonian H % per volume V 1 of the i-th cell can be 
expressed as a function of the z-th value of x l ; 

n.n. 

El,- !'-<■' i:.,, ■ /'-'•'• (i) 

The first self-energy term is quite similar to the lead- 
ing P 2 term in the Landau-Devonshire equation, which 
provides an approximation of the barrier height between 
the bi-stable points in the FE potential [2lJ . The higher 
order terms P 4 and P 6 were ignored here for the simplic- 
ity of the calculations. The second term describes the 
short-range wall interaction between the nearest neigh- 
boring domains, where r n . n , and x n , n . are the adjacent 
circumference and unit vector of nearest neighbor domain 
cells, respectively. It has a maximum value when the 4 
nearest neighbor domain cells surrounding a nucleated 
cell have the opposite polarization. The third term orig- 
inates from the long-range dipolar interaction between 
cells. The dipolar field E d is the sum of the dipolar elec- 
tric field caused by all the other cells. The last term 
describes the work done by external source. A similar 
calculation method was used in micro-magnetic systems 

E2 

When two electrodes sandwich a FE film, the electrode 
screens the value of E d produced by P s . The screening 
can occur not only in the perpendicular direction but 
also in the lateral direction. A numerical estimate re- 
veals that the lateral screening length should be of the 
order of the FE film thickness. For an ultrathin FE film, 
only the dipolar interaction between the nearest neighbor 
domains across domain boundaries could be important, 
so E d can be approximated to e~d J2n n r n.n. x n.n.' wnere 
id is the line density of the dipolar field across the cell 
boundary. Then, Eq. (1) can be rewritten as 

n.n. 

—E ex t ■ P s x l , (2) 

where an effective wall energy density &w can be de- 
fined as aw — cr w + edP s 1V % I dsTO- Note that the 
analytic form of the screened dipolar interaction energy 
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FIG. 4: (color online). A scaling relation between n and 
U* jksT. The open and solid symbols come from the experi- 
mental AP(t) data with and without E ex t, respectively. 



term is identical with that of the domain wall energy. 
The simulations were performed for room temperature 
domain dynamics with K 2 = 3xl0 6 J/m 3 , (Jw^= -1 
mJ/m 2 , and an attempt frequency of 10 9 Hz. 

The Monte-Carlo simulations were performed with and 
without lateral screening of the interaction term between 
long-range domains in the Hamiltonian. Without the lat- 
eral screening of the interaction term, the simulations for 
the 5.0 nm thick BaTiC>3 film predicted a very rapid ex- 
ponential decay of AP, as shown by the black dotted line 
in the inset of Fig. 3(a). When the screened domain in- 
teraction term was included as in Eq. (2), the simulations 
showed that the polarization decay would slow down sig- 
nificantly, exhibiting a power-law behavior, as shown in 
the inset of Fig. 3(a). The calculated polarization decay 
agrees with the experimental results quite well. Without 
changing any other parameter values, we performed the 
same simulations for other thicknesses of BaTi03 capac- 
itors, and found reasonably good agreement, as shown 
in Fig. 3(a). We also performed similar simulations for 
the cases with non-zero values of E ext , and obtained good 
agreement with experimental data, as shown in Fig. 3(b). 
These results clearly demonstrate that domain dynamics 
must be governed by thermally-activated nucleation pro- 
cesses in ultrathin films |2(ij . 

When the thermodynamic domain nucleation is the 
governing process for polarization decay, U* /ksT should 
be the most important physical term which determines 
the decay rate. Following the methods in Ref. ^3], we 
were able to convert the values of (\Ed\ - \E ext \) into 
U* values and made a plot of n vs. U*/kBT. Figure 4 
demonstrates that there might exist a universal relation 
between n and U^/ksT, regardless of dsro and (\Ed\ 
- \E ext \). To confirm the thermodynamic nucleation- 
governed polarization switching process, we indepen- 
dently performed polarization decay measurements for 
the 30 nm thick BaTi03 capacitor at various tempera- 



tures. The solid red circles in Fig. 4 correspond to the 
experimental data for temperature-dependent polariza- 
tion decay, which fall onto the seemingly universal line 
between n and U* /ksT. This indicates that there could 
be a simple relation between n and U* /ksT, which is uni- 
versal, irrespective of film thickness, total internal elec- 
tric field, and temperature. 

Our findings about the thermodynamic nucleation- 
governed switching process and the empirical relation be- 
tween n and U* /ksT have some important implications 
for the design of nanoscale FE devices. One implication 
concerns the capacitor geometry of FE devices where the 
polarization decay due to the Ed is predictable by the 
empirical relation we found. It could impose a new funda- 
mental thickness limit on such FE devices. For example, 
ferroelectric random access memories require the reten- 
tion of a certain value of AP for 10 years 0- For our 
BaTi03 capacitors, U* should be higher than 40 fc^T 
to retain AP at more than 50% of the initial value (n 
< 0.02). This value of U* corresponds to a thickness of 
about 40 nm for our fully-strained BaTi03 films, which is 
much thicker than the critical thickness (2.4 nm) for the 
ferroelectricity of BaTi03 [2j . For ultrathin film capaci- 
tors, a condition of U* > A- ksT, where A is a constant 
of the order of 10, has to be satisfied in order to have 
a stable retention property. As the Ed value increases 
(decrease in U*) with the reduction of (Ibto, retention 
loss due to the thermodynamic nucleation of reversed do- 
mains becomes much severe. Note that this situation 
is analogous to that of the superparamagnetic limit for 
magnetic memory devices: a long-range ferromagnetic 
order vanishes, when anisotropy energy becomes compa- 
rable to thermal fluctuation energy j23j . 

Our work has clearly shown the importance of a ther- 
modynamic domain nucleation process in the polariza- 
tion decay of ultrathin FE films. The polarization de- 
cay should impose a fundamental thickness limit on 
capacitor-type devices, which is much stricter than the 
critical thickness limit. Further experimental and the- 
oretical studies are clearly required to clarify the ori- 
gin of the empirically-determined relation between n and 
C/*/fcsT and its applicability to other ultrathin FE sys- 
tems. 

We thank M. W. Kim for discussion. We acknowledge 
financial support by the Korean Ministry of Education 
through the BK21 projects and the KOSEF through the 
CRI Program. 



* Electronic address: twnoh@snu.ac.kr 
[1] Y. S. Kim et al, Appl. Phys. Lett. 86, 102907 (2005). 
[2] J. Junquera and Ph. Ghosez, Nature 422, 506 (2003). 
[3] D. D. Fong et al., Science 304, 1650 (2004). 
[4] K. J. Choi et al, Science 306, 1005 (2004). 
[5] H. N. Lee et al., Nature 433, 395 (2005). 



[6] I. Kornev, H. Fu, and L. Bellaiche, Phys. Rev. Lett. 93, 
196104 (2004). 

[7] Bo-Kuai Lai et al, Phys. Rev. Lett. 96, 137602 (2006). 
[8] T. Tybell, P. Paruch, T. Giamarchi, and J.-M. Triscone, 

Phys. Rev. Lett. 89, 097601 (2002). 
[9] J. F. Scott, Ferroelectric Memories (Springer-Verlag, 

New York, 2000). 
[10] Y. W. So, D. J. Kim, T. W. Noh, J.-G. Yoon, and T. K. 

Song, Appl. Phys. Lett. 86, 092905 (2005). 
[11] R. Landauer, J. Appl. Phys. 28, 227 (1957). 
[12] A. Aharoni, Introduction to the Theory of Ferromag- 

netism (Clarendon Press, Oxford, 1996). 
[13] A. M. Bratkovsky and A. P. Levanyuk, Phys. Rev. Lett. 

85, 4614 (2000). 
[14] G. Gerra, A. K. Tagantsev, and N. Setter, Phys. Rev. 

Lett. 94, 107602 (2005). 
[15] D. J. Kim et al, Phys. Rev. Lett. 95, 237602 (2005). 
[16] Y. S. Kim et al, Appl. Phys. Lett. 88, 072909 (2006). 



■5 



[17] R. R. Mehta, B. D. Silverman, J. T. Jacobs, J. Appl. 
Phys. 44, 3379 (1973). 

[18] H. F. Kay and J. W. Dunn, Phil. Mag. 7, 2027 (1962). 

[19] W. J. Merz, Phys. Rev. 95, 690 (1954). 

[20] Our coercive field data agrees well with the prediction for 
the crossover from the half-prolate spheroidal to cylin- 
drical shapes at oIbto ~ 10 nm and also supports the 
picture of polarization switching dynamics governed by 
thermodynamic nucleation process in ultrathin ferroelec- 
tric films (to be published). 

[21] N. A. Pertsev, A. G. Zembilgotov, and A. K. Tagantsev, 
Phys. Rev. Lett. 80, 1988 (1998). 

[22] S.-B. Choe and S.-C. Shin, IEEE Trans. Mag. 36, 3167 
(2000). 

[23] S. Chikazumi, Physics of Ferromagnetism (Oxford Univ. 
Press, New York, 1997). 



